function [x,fval] = dispatch(A)
[m,n] = size(A);
c = A(:);
% 约束条件
a = zeros(m+n,m*n);
% 所有变量都是整数（只存在发配或不发配）
intcon = 1:m*n;
% 生成约束条件的系数
for i=1:m
    a(i,(i-1)*m+1:m*i) = 1;
    a(m+i,i:m:m*n) = 1;
end
%资源向量
b = ones(2*m,1);
lb = zeros(m*n,1);
up = ones(m*n:1);
[x,fval] = intlinprog(c,intcon,[],[],a,b,lb,up);
% 重新用矩阵显示x
x = reshape(x,[m,n]);
end
